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Abstract 

Periodic driving can create topological phases of matter absent in static systems. In terms 
of the displacement of the position expectation value of a time-evolving wavepacket in a closed 
system, a type of adiabatic dynamics in periodically driven systems is studied for general initial 
states possessing coherence between different Floquet bands. Under one symmetry assumption, 
the displacement of the wavepacket center over one adiabatic cycle is found to be comprised by two 
components independent of the time scale of the adiabatic cycle: a weighted integral of the Berry 
curvature summed over all Floquet bands, plus an interband coherence induced correction. The 
found correction is beyond a naive application of the quantum adiabatic theorem but survives in 
the adiabatic limit due to interband coherence. Our theoretical results are hence of general interest 
towards an improved understanding of the quantum adiabatic theorem. Our theory is checked using 
a periodically driven superlattice model with nontrivial topological phases. In addition to probing 
topological phase transitions, the adiabatic dynamics studied in this work is now also anticipated to 
be useful in manifesting coherence and decoherence effects in the representation of Floquet bands. 
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I. INTRODUCTION 


Just like the energy bands of an electron moving in a crystal, Floquet (quasi-energy) 
bands capture almost all aspects of a periodically driven system with translational symme¬ 
try. As shown by recent theoretical and experimental studies HH26], Floquet bands may 
possess many intriguing topological phases that are absent in static systems. It is impor¬ 
tant to connect Floquet band topology with observables and identify clear signatures of 
the associated topological phase transitions. To that end one fruitful topic is topologically 
protected edge states, but the issue of bulk-edge correspondence in driven systems is still 
under development jT2J EH :22, [233 • 

Thouless’s adiabatic pumping shares the same topological origin as the 2-dimensional 
integer quantum Hall effect [2ZH3Q]. There the sum of the Chern numbers of all filled en¬ 
ergy bands below the Fermi surface determines the number of pumped charges [27j (though 
not relevant to this work, we mention that Thouless’s seminal result is also applicable to 
many-body systems [25]). In non-equilibrium situations such as driven closed systems, there 
is no longer a Fermi surface to guarantee either filled or empty Floquet bands. Neverthe¬ 
less, an analogous link between Floquet band topology and adiabatic transport was also 
found in driven systems 0 E3]. Indeed, for a closed system (here we do not consider a 
system attached to leads as in a traditional transport problem), Thouless’s result may be 
re-interpreted as the quantized displacement of a band Wannier state over one adiabatic 
cycle [S], an interpretation fully consistent with the modern polarization theory [32] • hr 
parallel to this, the displacement of the position expectation value of a coherent wavepacket, 
prepared initially as a Wannier state of one Floquet band, can also be connected with the 
Floquet-band Chern number mm- However, this result is still of limited use in detecting 
topological phase transitions, because the initial-state preparation already needs the full 
knowledge of a Floquet band (which is often unavailable) and the initial state has yet to be 
re-prepared when a system parameter changes. 

This work considers a type of adiabatic dynamics in periodically driven systems for gen¬ 
eral initial states. To be as simple as possible our calculations are based on single-particle 
wavepacket dynamics, though the extension to many non-interacting particles initially oc¬ 
cupying many different Floquet states can be also done without difficulty. One of our 
motivations is to address more realistic situations, where the dynamics emanates from an 
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easy-to-prepare initial state regardless of other system parameters. The initial state then 
naturally occupies many Floquet bands and thus constitutes a coherent superposition across 
different Floquet bands. As shown below, such interband coherence (1BC) has a rather sur¬ 
prising contribution to the adiabatic dynamics. In short, under one symmetry condition, 
the displacement of the time-evolving wavepaeket is given by two components: a weighted 
integral of the Berry curvature summed over all Floquet bands, plus an IBC induced correc¬ 
tion not captured by a naive application of the quantum adiabatic theorem. The correction 
is independent of the duration of an adiabatic cycle. That is, no matter how slowly an 
adiabatic protocol is executed, the same correction emerges. We shall also briefly discuss 
the implications of the found correction. 

II. POPULATION CORRECTION IN ADIABATIC FOLLOWING 

We start with a one-dimensional driven Hamiltonian Hp(x,t), where x and t represents 
coordinate and time. The driving period is assumed to be r, he., Hp(x,t + r) = Hp(x, t ). At 
any t, the driving held maintains the translational symmetry, with Hp(x + a, t) = Hp(x,t), 
where a is the lattice constant. On top of the periodic driving, parameter f3 is tunable, thus 
allowing for an adiabatic protocol to be considered later. Without loss of generality (3 is 
assumed to be a periodic parameter with period 2n. Define U(f3 ) as the Floquet (one-period 
time evolution) operator and \ip n ,k(/3)) as Floquet eigenstates normalized over one unit cell. 
\^n,k(/3)) are Bloch states characterized by quasimomentum k, with eigenphases oj n ,k{^)- 
That is, 

u(mn,M) = ( 1 ) 

The collection of eigenphases u> n ^(/3) for k G [— Tr/a, n/a] and for /3 G [0, 27t] form the 
nth (extended) Floquet band. The Floquet bands are assumed to be gapped. The phase 
convention of the Floquet eigenstates is chosen to be ( , 0n,fc(/?)|d'0n,fe(/?)/d / d) = 0, the so-called 
parallel transport condition. Under this convention, | ip nt k(/3 = 2tt)) = e~ lln ’ k \ iJ} n ,k(/3 = 0)), 
where 7 ,^ is the Berry phase. 

Consider then an adiabatic protocol f3 = f3(s) of duration Tr, where s = t/(Tr ) is the 
scaled time, with /3(0) = 0, and f3( 1) = 27t. For example, f3(s) = 2i rs indicates a linear 
sweeping in (3. Because the Hamiltonian is periodic in f3 with a period 2n, such a protocol 
implements an adiabatic cycle. For convenience each small increment in (3 is introduced at 
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multiple (j) periods of driving, i.e., at tj = jr or Sj = j/T. We define the dynamical phase 


j 

) = ^ ' ^n.k [,d ('Sj' )] ; 
3 ‘=1 


( 2 ) 


namely, the accumulation of the instantaneous Floquet eigenphases co n ,k(/3 ) from s = 0 to 
s = Sj. For very large T, s J+1 — Sj = 1/T can be taken as the differential ds, which then 
gives an alternative, but integral form of the dynamical phase: f2 ni fc(s) = T / Q s cu ni fc[/9(s)]ds. 
As a somewhat standard technique, we decompose a time-evolving state l'h(s)) using the 
instantaneous Floquet eigenstates, 

l*M> = (3) 


where C n ^(0 ) depicts the initial state |(0)), with p n ,k(0 ) = |C ni fc(0)| J being the initial 
population on the nth band with quasimomentum k. Note also that each k component 
can be considered separately because k is conserved throughout. Now by our construction 
|tf(s + ds)) = U[/3(s + ds)]|T(s)). Projecting |\l/(s + ds)) [in connection with Eq. Q] onto 
("0 n,k[/3(s + ds)] | and also using our phase convention, one arrives at [33] 


dC n , fc 

ds 



^m,k) 


m,k 


m^n 



ds / * 


(4) 


Equation (|f]) describes possible transitions between instantaneous Floquet eigenstates of 
different band indices m and n. A naive application of the adiabatic theorem would be 
equivalent to setting dC/^/ds = 0. In the same spirit of Thouless’s derivation of adiabatic 
pumping [2T[ [28], we go beyond this crude adiabatic approximation by integrating Eq. Q 
and keep terms up to the order of 1/T. This results in (see [33] for details) 


Cn,fc(l) — C n ,k(0 ) + — 'y [ C m>k {0) | W nm ,k{s) 


m^n 


S— 1 ' 


s=0 , 


(5) 


where 


W n 


r,k( S ) = 


dpn,k(P) 

dlpm,k(/3) \ 

ds j 


[P)-Um,k(t 3)] 




( 6 ) 


Technically, Eqs. (J5|) and (|6| are not completely parallel to what can be expected from a 
conventional first-order adiabatic perturbation theory for non-driven systems. This is mainly 
because here, as s —> s + ds, the minimal change in the dynamical phase Q n:k {s) is given 
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by uj nt k[/3(s)\ [see [33] and Eq. ([2])]. Indeed, to our knowledge, Eqs. (|5j and § represent 
the first explicit result of a first-order adiabatic perturbation theory applied to periodically 
driven systems. 

Of more significance is the slight correction to the population on each instantaneous 
Floquet eigenstate, which cannot be identically zero in any finite-time adiabatic protocol. 
Let the population change on the nth instantaneous Floquet band with quasi-momentum k 
be 

Ap^k = \C n , k (l)\ 2 -p n ,k(0). (7) 

Then to the first order of 1/T, one obtains from Eq. (J5]) 
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j 

II 

A 

Zj 

^ Pn,k 

£ c;*(o)cv(o) w nm , k 


_rriy^n \ 

s=° / _ 


Two observations can now be made. If, as in previous work 0 123, the initial state is 
deliberately prepared on one single Floquet band, then the cross term C* k (0)C m ^(0) is 
zero and hence A p n ^ at least scales as 1/T 2 . This scaling is well known. By contrast 
and much less discussed, for a general initial state as a natural superposition state across 
multiple bands, the coherence cross term C* k (0)C mj k(0) is in general nonzero and as a 
result scales as 1/T. Such two types of scaling are caused entirely by IBC. That is, as 

far as the population correction is concerned, the adiabatic following of a superposition of 
instantaneous eigenstates is markedly different from that starting from a single instantaneous 
eigenstate. 


III. WAVEPACKET DYNAMICS DURING AN ADIABATIC CYCLE 


In this work the adiabatic following dynamics is considered in terms of the change in 
the position expectation value over one adiabatic cycle, denoted by A(x). Lengthy but 
straightforward calculations [33] show that, apart from a transient effect and for large T, 


A (*> = Z) / 

n 


d k 



d'T n,k 

“dF^ d k 


|C^(1)| 2 , 


d r /n,k . 

d k 


[^Pn,k d" Pn,/c(0)] 7 


(9) 


where the role of the above-discussed population correction A p n ^ is highlighted. Upon a 
multiplication of the two factors in the integrand of Eq. g. four terms emerge. Because the 
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geometric phase 7 is independent of T and A p n ^ scales as 1/T, the first term p H)k 

vanishes as T —>■ + 00 . The second term ^^-p n ,k( 0) is the Berry phase derivative weighted by 
an initial distribution. Of special interest is the third term dfln d ^ Ap n ,k- Even though A p n ^ 
scales as 1/T, the dynamical phase 0^(1) and hence its derivative is proportional 

to T. As such this product scales as 1/T° and will not vanish even when T —y + 00 ! Finally, 
the fourth term which reflects the influence of ballistic motion on A(x), for 

an initial distribution p n ,k{ 0)- With some symmetry in the spectrum and in p n ,k( 0), this 
fourth term vanishes upon integration of k. For example, a ^-reflection symmetry in the 
populations p n ,k{ 0) = p n ,-k( 0) and in the Floquet spectrum cu n _ k (/3) = co nt k(/3) suffices. 
Specifically, such kind of symmetry is often the case if the initial state and the nondriven 
version of the system has left-right symmetry. 

In the adiabatic limit and under the above symmetry assumption, Eq. reduces to 


A ( x ) = 5Z / 


d k 


d 7n,fc ^ , df2 nj fc(l) A ^ 

Pn,k\ 0) T .. Ap nk 


d k 


d k 


( 10 ) 


It is useful to rewrite as f d/3B n (/3, k), where B n (/3,k ) is the Berry curvature of the 
nth band, with 


B n {P,k) 


/dlJ> n ,k{P) |^n,fc(^)v /<%»,*(£) I &tp n> k(P)\ 
' dk 1 d/3 ' ' <9/3 1 dk ' 


( 11 ) 


It is also time to use the expression of A p n ^ from Eq. (J8]) . Upon integration over k, only the 
W nro , fc (0) part in A p U)k survives. To conclude our theory, we dehne the average quasienergy 
E Ujk along a pumping cycle, i.e., 


E n , k = (l/T)Q nt k(l) 


w n , fc [/3(s)]ds. 


( 12 ) 


One then arrives at 


A(x) = ^2 I dk J d/3 B n (/3,k)p nt k( 0 ) 


m^n * 


dk Re 


c;yo)c m vo)^MW(o) 


(13) 


Both components in Eq. (13) are independent of T. One is just a Berry curvature integral 
over a 2-torus of (/3, k) with a nonuniform weighting factor p n ,k( 0) (A somewhat similar for¬ 
mula was given in Ref. [3] regarding conductance under ac and dc fields). The other compo¬ 
nent is the IBC induced correction. I 11 addition to “off-diagonal” coherence C* k (0)C mi k(0), 
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FIG. 1: (color online) (a) Floquet eigenphases of CDHM plotted as a function of k and (3. (b) 
Probability distribution function p n *.( 0) on three bands, for an initial state exclusively located at 
site l = 0. The CDHM parameters are J = K = 3. 

the correction is also related to lFnm,fc(0). As indicated from Eq. (J6|, is proportional 

to 1/ {l — e*t a;ri ’ fc ^ 0 ' )_a ' m ’ fc ( 0 )]}. So if two Floquet bands are almost touching at a particular 
value of k, then this correction itself may become singular. Therefore, the correction itself 
may be useful in detecting a topological phase transition. Furthermore, because PFnm^O) 
presented in Eq. (J6]) is proportional to d/3(s)/ds at s = 0, it is seen that the found correction 
depends sensitively upon how the adiabatic protocol is turned on in the very beginning. 
In this work the switching-on of the adiabatic protocol is characterized by d/3(s)/ds = 2n, 
but our theory applies to other situations as well. It is also interesting to note that, if the 
initial state is made to uniformly occupy one single Floquet band (hence no IBC), then the 
correction term is zero. At the same time the uniform Berry curvature integral precisely 
gives the Floquet band Chern number. In retrospect, previous treatment mm based on an 
idealized adiabatic theorem works simply because, without IBC the population correction 
A p n k scales as 1/T 2 instead and therefore it can be safely neglected. As a final remark, 
we stress that if the initial state is not fixed and injected randomly ( e.g in a transport 
experiment where the system is connected with two leads), then the IBC term needs not to 
be considered because it will self-average to zero. For this reason we emphasize here that 
our results are for a closed system that is periodically driven. 
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IV. THEORY VS MODEL CALCULATIONS 


We apply our theory to a continuously driven Harper model (CDHM) [i23 t [23]. CDHM 
describes a particle hopping on a periodically modulated superlattice with its dimensionless 
Hamiltonian given by 


H = ^(J/2)(a ^ +1 + h.c.) 

i 

+ K cos(27 rt/r) cos (2nal + j3) aja/, 


(14) 


where / is the lattice index, aj (a/) is the creation (annihilation) operator, J is the hopping 
amplitude, and K is the potential strength. This model may be potentially realized by 
optical waveguide experiments or by cold atoms in a deep optical lattice [21, 25], In accord 
with our general notation above, the adiabatic parameter is /3, with /3(s) = 2irs. For 
a = M/N (M and N are two co-prime integers), the superlattice potential has a period of 
N, yielding N well-gapped Floquet bands in general. Here we choose M — 1 and N = 3, 
with r = 2. The initial state is fixed, which is placed exclusively at the site l = 0. Imagine 
an optical waveguide realization of our model. Then onr initial state requires that initially 
the light is injected in one waveguide located at site zero. For a cold-atom setup, our initial 
state can describe an ensemble of bosonic atoms initially all at site zero. Figure [l](a) shows a 
typical Floquet band structure in the case of J = K = 3. When analyzed in terms of Floquet 
bands, our seemingly simple initial state populates all three bands. Figure [ljb) depicts the 
respective initial population on the three Floquet bands, which is denoted as p n ,k( 0) in our 
theory. Note also that both the spectrum and the initial state satisfy a reflection symmetry 
in the k space. 

Presented in Fig. 2 is the population change A p n ^ vs k after one pumping cycle, for 
the bottom Floquet band (see Fig. 1). Panel (a) displays the actual A p n ,k obtained from 
dynamics calculations. Panel (b) depicts our theoretical result from Eq. ([8]) based solely on 
an IBC analysis. The agreement between theory and numerical experiment confirms our 
adiabatic perturbation theory for periodically driven systems. 


Next we check our central theoretical prediction presented in Eq. (13). To that end we 
consider the same adiabatic protocol with a varying duration T. As shown in Fig. [3j the 
one-cycle displacement in the wavepacket center stays the same as we change T over a wide 
range. The actual result is also compared with the integrated Berry curvature weighted by 




FIG. 2: Actual (a) and theoretical (b) population change A p n ^ vs k, after one pumping cycle in 
CDHM, for J = K = 3, and T = 1024. 


p n ,k (0) (triangles in Fig. [3]). It is seen that they are far away from each other. Finally, 
the IBC effect beyond a naive application of the quantum adiabatic theorem [presented in 
Eq. (13)] is included, with the corrected prediction represented by the filled circles in Fig. [3j 
The theory perfectly agrees with our dynamics calculations. 

We now investigate the adiabatic dynamics in the vicinity of a topological phase transi¬ 
tion. We record and present A(x) in Fig. 4 while scanning the value of J = K 6 [5.0, 5.3]. 
The actual result, smoothly varying over almost the whole shown regime, becomes highly 
erratic in a small window around J — K — 5.15. This behavior is fully consistent with 
the jump of the Floquet-band Chern numbers, from (4,—8,4) to (—8,16,—8) at about 
J — K — 5.14 (earlier we observed this topological phase transition in the same model in 
Ref. [23]). Our theoretical result also displays a sudden jump at about J — K — 5.14. 
In the immediate neighborhood of this sudden jump, the actual result becomes sensitive 
to T (not shown) and deviates from our theory based on adiabatic following. This behav¬ 
ior is expected because there the Floquet bands are almost touching each other, yielding 
more nonadiabatic effects for a protocol of a finite duration. We stress that the successful 
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FIG. 3: Time dependence of the position expectation value during one adiabatic cycle in 
CDHM. The six solid lines (ending with empty circles) describe the actual dynamics result, for 
T = 1024,2048,3072,4096,5120,6144. J = K = 4. Dashed line connecting filled circles represents 


our theory in Eq. (13). Dashed line connecting triangles represents the Berry curvature component 
only. 


detection of a topological phase transition illustrated here is achieved using a very simple 
initial state regardless of system parameters. Note also that the sole term of the weighted 
Berry curvature integral (triangles in Fig. 4) is far away from the actual displacement of the 
wavepacket center. 


V. SUMMARY 

Under one symmetry assumption, a type of adiabatic wavepacket dynamics in a closed 
and periodically driven system is studied. The associated displacement of wavepacket center 
over one adiabatic cycle is found to be comprised by two components independent of the 
duration of an adiabatic protocol. In addition to a weighted integral of the Berry curvature, 
there is an IBC induced correction. This correction survives in the adiabatic limit, making it 
an important ingredient in an improved understanding of the quantum adiabatic theorem. 
The explicit expression of the found correction also indicates that it depends sensitively 
upon how fast an adiabatic protocol is switched on in the very beginning. If we regard 
the adiabatic protocol as a result of a very slow driving held, then the observation that 
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FIG. 4: Detection of a topological phase transition in CDHM with J = K. The actual wavepacket 
displacement (empty circles, for T = 2560) agrees with our theory (filled circles) except for a 
window around J = 5.15. Triangles are the results excluding the IBC induced correction. 

how an adiabatic protocol is switched on can affect the dynamics over an adiabatic cycle is 
somewhat consistent with the following: in driven systems the switching-on protocol of a 
driving field does play an important role [35]. 

Because the wavepacket displacement studied here depends on quantum coherence in the 
representation of Floquet bands, our theory is also of interest to coherent control studies [33] 
that advocated the utilization of quantum coherence to influence quantum dynamics. Indeed, 
one may further manipulate the IBC effect (e.g., by considering adiabatic protocols different 
from what is studied here) and thus achieve a coherence-based control of the adiabatic 
dynamics. In addition to probing topological phase transitions, the adiabatic dynamics 
studied in this work is anticipated to be useful in manifesting coherence and decoherence 
effects in driven systems. To that end, a recent study of dissipative Floquet topological 
systems [35] can be a useful starting point. 

Acknowledgments: We thank Derek Y. H. Ho for helpful discussions. 


In this Appendix, we use the same notation as in our main text. The first part is about a 
detailed derivation of our first-order adiabatic perturbation theory (he., to the order of 1/T) 
for periodically driven systems under an adiabatic pumping protocol. In the second part 
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we present a detailed derivation of Eq. (8) of the main text, namely, how the change in the 
position expectation value of the system over one adiabatic cycle is related to the dynamical 
phases, geometric phases, as well as the final state population distribution introduced in the 
main text. 


Appendix A: First-order adiabatic perturbation theory in periodically driven sys¬ 
tems 


As stated in our main text, we assume that each slight change in the adiabatic parameter 
f3 is introduced at multiples of r, where r is the driving period. This assumption is not 
essential, but it brings us much convenience when we describe a time-evolving state as a 
superposition of instantaneous Floquet eigenstates, denoted as (/?)). Consider then an 
adiabatic protocol of duration Tr. The protocol is summarized by f3 = f3(s), where s is the 
scaled time, with /3(0) = 0, and j3( 1) = 2tt. In particular, after j driving periods, t = jr, 
s = Sj = j/T. For very large T, Sj +1 — Sj = l/T can be regarded as the differential ds. 

Next we expand the evolving state |T(s)) as the following: 

I - 7T 

!*(.» = t/A /“ , (Al) 

_ a Tl 

where (s) is the accumulated dynamic phase, C^^O) depicts a general initial condi¬ 
tion |T(0)), with p n ,k( 0) = |C' n; fc(0)| 2 being the initial population on the n-th band with 
quasimomentum k. According to this construction, we have 


|T(s + ds)) = U\/3(s + ds)]|T(s)), 


(A2) 


where £/[/3(s + ds)] is the one period evolution operator defined in the main text. Projecting 
|T(s T ds)) onto the bra state (t/yfc [/3(s + ds)] | and using our parallel transport phase 
convention defined in the main text, one immediately has 


C n , k (s + ds)e- iQ ^ s+d ^ = e --n,.(s+ d S ) J2Cm,k(s)e- inm ’ h{s) (i>n,k [P(s + ds)] \^ m , k [/3(s)]> . 

m 

(A3) 

Note that k is conserved so different k components are not coupled to each other. The right 


hand side of Eq. (A3) can be divided into a term with m = n and other terms with m ^ n. 
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Moving the m — n term to the left hand side of Eq. (A3), one finds 


^c n , k (s) = - [£(s)] \^ m ,k [/3(s)]> • (A4) 

m^n 

To proceed we take special note of the actual meaning of d [ e _ *W,fcb)]. It means the 
change in the dynamical phase factor after one more slight change in (3 is introduced upon 
one more period of driving. This indicates that the minimal change in the dynamical phase 
as s —> s + ds is uj Uy k[(3{s + ds)], i.e., the instantaneous eigenphase of the Floquet operator 
at f3(s + ds). Therefore, 


d g iQn,k (s+ds) g __ [/3(sH-ds)] g 


This is to say, 

e _ * n ™’ fc ( s )ds = 

Via the same reasoning, we have 

J [-Ti, A- ('5') — 


l/T 


e -*^n,fc[| 3 (s+ds)] _ 

l/T 


d [e 


(A6) 


*{ aj ".'=[/ 3 ( s + ds )] _aj m.'=[/ 3 ( s + ds )]} _ 


(AT) 


We are now ready to integrate Eq. (A4) by parts via the relation in Eq. (A7). We hence 
obtain 


C n ,k( s ) = 


C, JO) 4- — Cm,k(s )(j>n,k [ff(s )] | ds ;^m,fc [ff(s')]) 


T 

1 

T 


m^n 


s '=0 


/ f ] 6 ./ ^[W.fc^O-W.fcbO] d / Cm,k(s ){jj n ,k [/3(g )] \ ds ,j>m,k [/3(g )]) 
J I- !«' ' X — e H Wn ' fe [ , 8 ( S, )] _a, m 1 fc[l 8 ( s ')]} 


m^n 


(AS) 


The third term on the right hand side of Eq. (A8) is at least of order l/T 2 , that is: 


1 

J 1 '! 


T ^ 

rriy^n ’ 


E 

m^n _ 


4 S ' e *[ n n,fc(s')-n m , fc ( S ')]_^_ ( Cm,k(s )(Vjfc [/3(s )] | dg/VVfc [/3(s )]) 
ds' l l _ e *{w„,fc[y9(a , )]-w m ,fc[j8(s')]} 

e i[n„,J S ')-a m ,Ja')] d {C m ,k{J)($n,k\P(s')]\T?'tl>m,k\fi(8')]) 


l — e '{ w r i ,#(s')]- u m,t[/ 3 ( s ')]} ds' l _ e *{^n,fc[^(s')]-^m,fe[^(s')]} 




m^n ' 


d ( Cm,fc(s , )(W,fe[/ 3 (s')]l)fJ 7 V’m,fe[/ 3 b , )]> 


ds' y l e *{“n,^( s ')]-“m,fe[W s ')]} 


s'=s 

s '=0 


(A9) 
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Ignoring this term of the order of 1/T 2 , we find that up to order 1/T, the amplitude on the 
nth band with quasimomentum k at scaled time s is given by: 


c n , k (s) = c nr , fc (o) + ^ ^ c m , k ( o) [w nm , k (s')} i;!:; 


(A10) 


m^n 


with 


) — 


(All) 


X — [/S(s')] —^TTJ.,fc [/3(-5 , )l } 

The above two equations are just the results reflected by Eqs. (5) and (6) in our main text. 
Note that here we still retain those terms of the order of 1/T because, as seen from our 
theoretical analysis, a combination of the population correction of the order of 1/T and 
a dynamical phase factor proportional to T will affect the adiabatic dynamics even in the 
adiabatic limit. 


Appendix B: Adiabatic dynamics: Change in position expectation value 


Here we present a detailed derivation of Eq. (8) of the main text. The goal is to connect 
the change in the expectation value of position with a weighted integral of the Berry phase 
derivative and the dynamical phase derivative flUr V/ ( 1 ' . The weighting factor |C ni fc(l)| J 
in Eq. (8) of the main text is the final state population on the nth band with quasimomentum 
k. Throughout k G [— -, -] is assumed to be in the first Brillouin zone, where a is the lattice 
constant of the driven system. 

Qualitatively, Eq. (8) of the main text is somewhat intuitive because the position operator 
x should be somewhat related to i-A However, for a general initial state involving interband 
coherence (IBC), it is helpful and useful to present a detailed derivation. This task is 
straightforward but also quite technical. We start by rewriting the Floquet Bloch states 
| ip n! k) in terms of \u n ^) normalized over one unit cell, with u n ^(x) = (x|u n ,fe) being a 
periodic function of x with period a. That is, 


(Bl) 


(a#n,k) = e lkx (x\u n , k )] 

(X T (l\Un,k) (p^\'U"n,k) • 

Consider now an arbitrary state T), expanded as a superposition state across different 
Floquet bands, i.e.,, 

l») = f_[dkJ2 D n,Mn,t) ■ (B2) 
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Unlike Eq. (1) in the previous section or Eq. (3) of the main text, here we do not separate out 
a dynamical phase factor from the expansion coefficient. The (3 dependence of the Floquet 
states is also suppressed here because we first target at an expression for an arbitrary value 
of f3. The position expectation value on this state is given by ( x) = (T|£|T). In particular, 


(x) = 7jT dec d k d k! ^ DnA^n^x) X An,fc'(3#m,fc') 


2 7T 

a 

27T 

a 

2tt 


/ oo /*— n — 

dx J a dk J a dk' D r, 


k u* nk {x)e l " x D mk ,u mk i (x) x e 


ik x 


(B3) 


dx / ik ££>;,*<, 


—ikx 


dk' D ni}k 'U m:k '(x) ( -i —7 I e lk ' x 


d 


Performing an integration by parts in the last line of the above equation and using the fact 
fe'=s 


[D mtk iu m ^(x)e lk ' x ] 


k'=- n - 


= 0, one immediately has 


dk' D nhk iu m , k '(x) ( ) e zk ' x = i / dk' e lk ' x -^-D m yu, mik '{x) . (B4) 


d 

dk' 


Substituting the relation of Eq. (B4) back into Eq. (B3), we then obtain 




J[dkJ[dk' £ K,k J~ d* e-‘ (k - k> ul h (x)-^D m: k ,u m Ax) ■ (B5) 


To proceed further, we separate the integration over x into a summation over a unit cell 
index j and an integration over one unit cell. That is, 






-i(k-k’)(x+ja) * 

e a n,k 


{x) 


d_ 
dk ' 


Dm^'lLm^k’ (*^) 


(B6) 


r a h — k' - r) 

= dx 5 ( 2tT a + 3) e ^ (k ~ k ' )X KA^Q^ D m,k' u m,k'(x) , 

0 j =~ 00 

where x = x + j ■ a and x denotes the coordinate within one unit cell. In reaching the above 
expression the Poisson summation formula: 


OO OO j j j 

£ £ S(tx± a +j) (B7) 

j =—00 j =—OO 
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is also used. Now plugging Eq. (B6) into Eq. (B5), we find 


d 


{x) = i dx I dk ^Dl k u* njk (x)—D m , k u m , k (x) 

0 ~a n,m 

pa f— _ r 

die / dk ^ 


= i 


= i 


~ * d 9 

U n,k(^') U m,k{x)D n k~Qj^Dm,k + ^n,k^rn,kU n k {x^-^U m ^k{x) 


dk 


dk 


d 


d 


d 


'y 1 Dn,k cjj,Dn,k + ^ ^ |-Pn.,fc| ('Wra,fc | | w n,fc) + ^ ] ^n,k^fn,k ( u n,k \ ^_\ u m,k) 


dk 


n^=m 


dk 


(B8) 


Let us first focus on the third term in Eq. (B8). Because D n and D m ^ are the quantum 


amplitudes on different Floquet bands with different quasienergies, in general this term will, 
as time evolves, become a highly oscillating function of k and will then make a vanishing 
contribution to (x). To best illustrate this, it suffices to consider a driven system in the 
absence of the adiabatic protocol, be., at a fixed (3, with an evolution time Tr. In that 


simplified case, D n ^(s) = -D n ,fc(0)e * a V'=(£) T . Then the third term in the last line of Eq. (B8) 
reads as: 

/ — 

“ dk e~ iT ^ m -^ m ^G nm , k , (B9) 

a 

where G nni)k = D^ k (0)D mtk (0)(u nt k\-§ k u mtk )- It is now clear that this term will involve an 
integral of a highly oscillating function due to the factor e ~‘ lT [ ulrn ^^ Wn ^G)) j n the integrand. 
Analogous to the technique used in our adiabatic perturbation theory detailed in the previous 


section, one can now perform an integration by parts in Eq. (B9). It is then evident that 


this term is at most of the order of 1/T and can thus be neglected. Note that this third 
term does not have to be negligible at the start time due to the lack of the highly oscillating 
phase factor. Nevertheless, even if at the very start of an adiabatic protocol it happens to 
be nonzero, it is still a transient effect, insofar as at the start of the second adiabatic cycle, 
this term still gives a negligible contribution to (x). In addition, in our numerical studies of 
the continuously driven Harper model, this term is found to be already vanishingly small at 
the start of the first adiabatic cycle. 


After dropping the third term in the last line of Eq. (B8), let us now compare the difference 


in ( x) between the initial state and the final state of a pumping protocol. For the initial 
state defined in the main text and in the representation of |= 0)), we have 


Ai,fc — G nk (0). 


(BIO) 
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From Eq. (B8), one finds the initial position expectation value 


{x)(s = 0) = i / d k 


y ' ^rt.fc(O) ~XT^ , n,fc(0) + |Cn,fc(0)| {^n,k \ ~WT |^n,fc) 


(BH) 


where |u nj fc) is now associated with Floquet eigenstates for j3 = 0. For the final state, 
according to Eq. (3) of the main text, 



l*(s = 1)) = I ‘ d k J2C n AVe~ iQnM1) WnAP = !))■ 


(B12) 


Using | ijjn,k{P = 1)) = e lJn ’ k \i^ n ,k(P = 0)), the expansion coefficients -D n ,fc, still in the 
representation of | ^ n ,k{^ = 0)), is found to be 


D„, t = 


(B13) 


Using Eq. (B8) again, one finds the corresponding (final) position expectation value 

<*>(»= l) = i|‘di _ 

a \_ n n 

— r - ^dq n ,k dh2 ni fc(l)^ r-\\\2 


cq( 1) (1) + 5] I c„, t (1) I 2 (u„, fc I A |„ ni ) 


E 


d k 


n a 


+ 


(B14) 


d k d k 

The final result in terms of A(x) = {x)(s = 1) — {x)(s = 0) can then be obtained by 


subtracting Eq. (Bll) from Eq. (B14). As indicated by Eq. (A10) from the last section of 


this Supplementary Material, (7^(1) and C n ^{ 0) differ only by a term of the order of 1/T. 
Therefore, for sufficiently large T, this subtraction safely cancels most of the terms and leads 
to 


A (^) = 5Z 


+ ] |c n , t(1) | 2 . 


d k 


d k 


(B15) 


This is precisely Eq. (8) of the main text. Note that in Eq. (B15) we cannot replace |C' ni / c (l)| J 
by |<7 ni fc(0)| 2 as in treating other terms simply because it is accompanied by the factor 
that is proportional to T. 
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